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^ : ABSTRACT 
> ■ 

^ ' The negUgible eccentricity of all extra solar planets with periods less than 

^ ■ six days can be accounted for by dissipation of tidal disturbances within their 

00 ' envelopes which are induced by their host stars. In the period range of 7-21 days, 

^ ' planets with circular orbits coexist with planets with eccentric orbits. These will 

O ■ be referred to as the borderline planets. We propose that this discrepancy can be 

• attributed to the variation in spin-down rates of young stars. In particular, prior 
to spin-down, dissipation of a planet's tidal disturbance within the envelope of a 

2 ' sufficiently rapidly spinning star can excite eccentricity growth, and for a more 

^ ■ slowly spinning star, at least reduce the eccentricity damping rate. In contrast, 

• tidal dissipation within the envelope of a slowly spinning low-mass mature star 

■ can enhance the eccentricity damping process. Based on these results, we suggest 

■ that short-period planets around relatively young stars may have a much larger 
■ — ' dispersion in eccentricity than those around mature stars. We also suggest that 

because the rate of angular momentum loss from G and K dwarfs via stellar winds 
is much faster than the tidal transfer of angular momentum between themselves 
and their very-short (3-4 days) period planets, they cannot establish a dynamical 
configuration in which the stellar and planetary spins are approximately parallel 
and synchronous with the orbital frequency. In principle, however, such config- 
urations may be established for planets (around G and K dwarfs) with orbital 
periods of up to several weeks. In contrast to G and K dwarfs, the angular mo- 
mentum loss due to stellar winds is much weaker in F dwarfs. It is therefore 
possible for synchronized short-period planets to exist around such stars. The 
planet around Tau Boo is one such example. 
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1. Introduction 

Radial velocity surveys around nearby solar-type stars have led to the discovery of 
planets with short-period and highly eccentric orbits (Marcy et al. 2000). The period 
distribution is essentially logarithmic ranging from a few days to many years. While the 
upper range of the period distribution at this stage is incomplete due to relatively short 
survey times, there is a minimum cutoff at around period P — 3 days. The observed 
eccentricities, e, of planets with P > 21 days are uniformly distributed in the range zero to 
about 0.7. However, all planets with P < 6 days have nearly circular orbits (Figure 1(a)). 
In this paper, we consider the effect of star-planet tidal interactions on the eccentricity 
distribution of the short-period planets (hot Jupiters) . 

The dichotomy in the planets' eccentricity distribution is similar to that of binary stars 
(Mathieu 1994). It is generally assumed that the orbital eccentricities of binary stars with 
periods less than some critical value (the circularization period Pc) are damped by their tidal 
interaction (Zahn 1977). The magnitude of P^ is a function of the tidal dissipation efficiency 
and the age of the binary. The former quantity is usually parameterized in terms of the stars' 
(5*-values (cf Murray & Dermott 1999; in what follows, an asterisk subscript will refer to 
stellar values while the subscript p will refer to planets). Based on the observed magnitude 
of Pc in various stellar clusters (Mathieu 1994), the inferred values of Q'^ — 3Q/2k^ (where 
is the tidal Love number or twice the apsidal motion constant) are ~ 1.5 x 10^ (with 
Pc ^ 4: days) for young stars with age r* < 10^ yr, and ~ 10^ for mature stars (Terquem et 
al. 1998). 

The circularization of short-period planetary orbits may also be due to tidal dissipation 
processes (Rasio et al. 1996, Marcy et al. 1997). For example, this would be necessary to 
account for the small eccentricities of short-period planets if they were scattered into the 
vicinity of their host stars by other planets (Rasio & Ford 1996, Murray et al. 1998). In con- 
trast, if short-period planets acquired their small semi-major axes through tidal interaction 
with their nascent disks (Lin et al. 1996), planet-disk resonant interactions would naturally 
damp the eccentricities of planets with masses Mp < lOMj (where Mj is the mass of Jupiter) 
(Goldreich & Tremaine 1980, Artymowicz 1992, Papaloizou et al. 2001) although in some 
circumstances, it may also excite planets' eccentricities (Goldreich & Sari 2003). The pro- 
cesses (planet-star tidal interaction, unipolar induction, and disk depletion) which halted the 
orbital migration of hot Jupiters (Lin et al. 1996) may also induce the excitation of orbital 
eccentricity. In either early dynamical evolution scenarios, efficient post-formation eccen- 
tricity damping is needed to account for the small eccentricities (< 0.1) of the short-period 
planets. 

The orbital semi-major axes of short-period planets are larger than the radii of their 
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host stars by an order of magnitude and of their own radii by two orders of magnitude. The 
angular momenta of these planetary orbits are comparable to that carried by the spin of 
their host stars and larger than that contained in their own spins. Unless their orbits are 
synchronized and circularized, these planets are subject to intense tidal disturbances from 
their host stars. The Qp-vaAue for Jupiter has been inferred from lo's orbital evolution such 
that 5 X 10^ < < 2 X 10^ (Yoder and Peale 1981). With this range for Q'^, planets 
with periods less than 6 days whose masses and radii are comparable to those of Jupiter are 
circularized during the main sequence lifetime (r* ~ a few Gyr) of solar-type dwarf stars. 
This estimate is consistent with the observed circular orbits of the shortest period planets 
(TriUing 2000). 

The dissipation of these disturbances leads not only to dynamical changes but also 
to internal heating which may cause them to inflate (Bodenheimer et al. 2001). With a 
sufficiently large eccentricity, hot Jupiters would undergo runaway tidal inflation and overfill 
their Roche lobes before their orbits became circularized (Gu et al. 2003, 2004). Mass loss 
through Roche lobe overflow may either totally disrupt the planet or induce the orbital semi- 
major axis, a, to expand (Trilling et al. 1998). Both outcomes inhibit protoplanets from 
becoming too close to their host stars. Thus, the excitation and damping of eccentricity 
also has observable consequences and important implications for the structure and survival 
of individual planets and the period distribution of planetary systems. 

The observed eccentricity distribution also indicates that planets with periods between 
1 and 3 weeks have non-negligible eccentricities which arc smaller than those of planets with 
longer periods and larger than those of planets with shorter periods. We will refer to planets 
in this period range as borderline planets. In this paper, we adopt a conventional approach 
using frequency independent Q^^-values to analyze the orbital evolution of a star-planet 
system (Mardhng & Lin 2002). 

In §2, we show that the obscrvationally inferred eccentricities of these borderline planets 
are not well correlated with the expected time scale, Tgp, for eccentricity damping due to the 
dissipation of the stellar tidal disturbance inside the planets. Although it is customary to 
assume that dissipation inside a planet provides the dominant contribution to its orbital 
eccentricity evolution, dissipation within its host star can have a comparable effect if Q'^ ~ 
Q'p. We examine the possibility that, in this limit, the dispersion in the (e— Tep) relationship is 
caused by the observed spread in stellar rotation rates. Stellar spin evolves as a consequence 
of structural changes within a star, as well as mass loss and tidal interaction with a planet. 

In §3, we consider how these effects combined affect the long term evolution of a star's 
spin as well as the eccentricity of a planetary orbit. We show that short-period planets 
rapidly evolve into a spin-orbit quasi-equilibrium which continues to evolve, albeit at a 
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slower pace, because the dissipation of the stellar tidal disturbance within a planet continues 
to induce eccentricity damping. In contrast, during the early stages of its main sequence 
evolution, a solar- type star rotates rapidly and is likely to excite the eccentricity of its planet. 
However, as it loses angular momentum and attains a modest or slow spin frequency, the 
tidal dissipation within such a star will induce eccentricity damping. We also deduce the 
condition for quasi-equilibrium in which the loss of angular momentum through a stellar 
wind is balanced by that transferred to a star through its tidal interaction with a short- 
period planet. We then explore the possibility of such a system evolving into a state in 
which the stellar and planetary spins are aligned, and the planet and the envelope of the 
star spin synchronously with the orbit. Finally in §4 we summarize our results and discuss 
their implications for the dynamical evolution of short-period planets and the spin-down 
history of their host stars. 

2. The Eccentric Short-Period Planets 

Observational determinations of planetary eccentricities require extremely accurate ra- 
dial velocity data over many orbits. Spurious determinations of periods and eccentricities 
may be caused by sparse sampling, modest velocity measurement errors, and perturbations 
by other planets. Neglecting these potential sources of error in the published data, we find 
that half of the known systems with 1 week < P < 3 weeks have non negligible eccentricities. 
(We consider the present threshold of real detection to be e = 0.1 which corresponds to the 
typical ratio of residual to amplitude of radial velocity curves.) In this section, we consider 
the implications of these measurements. 

2.1. Dispersion in the Eccentricity Evolution of Short-Period Planets 

Using the standard formula (Goldreich & Soter 1966, Yodcr & Peale 1981, Peale 1999, 
Murray & Dermott 1999, for the limitation and uncertainties see §1), the timescale associated 
with eccentricity damping due to dissipation of the stellar tidal disturbance in a planet can 
be expressed as 

_ _ e _ ^ /^y g; / g; \ (m^\ /ma fj^\ [ra ~' 

(1) 

where is the mass of the star, a is the semi-major axis, and n is the orbital mean motion. 
The derivation of this formula assumes e <^ 1 as well as aligned and synchronous planetary 
rotation, that is, Qp — n — {GM^/a^y^'^. We will use it as a reference timescale to illustrate 
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the effects of significant eccentricity and stellar spin. 

In Figure 1(b) we plot eccentricity as a function of Tep of all known planets with P < 2 
months. The period range in this sample is selected to avoid the observational incompleteness 
for planets with longer periods. The magnitude of Q'p is assumed to be 10^ (see discussions 
in §2.3). The results in Figure 1(b) clearly show that e is negligibly small (< 0.1) for 
planets with r^p < 1 Gyr. In the transition range where 1 Gyr < r^p < 10 Gyr, planets on 
circular orbits (e < 0.1) coexist with those on eccentric orbits (e > 0.1). This dispersion 
cannot simply be attributed to the difference in the stellar ages (r^,) since the target stars 
are chosen for their lack of chromosphcric activity which generally indicates ~ a few Gyr 
(see Figure 1(c) for e versus Tep/r* ). In addition, their dynamical properties cannot be the 
main cause because some planets have quite different eccentricities despite the similarity 
between their periods and mass ratios. Although modest variations in Rp may significantly 
modify the magnitude of Tep, planets with P > 3 — 4 days cannot be easily inflated by tidally 
dissipation before their eccentricities are damped out (Gu et al. 2003). The insolent effect 
of stellar irradiation in quenching a planet's Kelvin- Helmholtz contraction (Burrows et al. 
2000) is also not well correlated with eccentricity (see Figure 1(d) for e versus stellar flux at 
a). If the short-period planets were scattered to their present location through dynamical 
instabilities on a variety of timescales, large eccentricities (~ 1) would be expected for all 
planets with Tep > and some planets with Tep smaller than r*. This inference is not 
consistent with the results shown in Figure 1(c). We cannot yet rule out the possibility that 
the dispersion in the e — Tep relation is due to some large- amplitude fluctuations in the Q'^ 
value which itself may be determined by some stochastic processes (see §2.3). 



2.2. Dissipation in the Host Star of the Planetary Tidal Disturbance 



The expression in eq(l) does not include the contribution from the dissipation in the 
star of the planet's tidal disturbance. In the case that both the stellar and planetary spins 
are ahgned with the orbit, the rate of change of e (Eggleton et al. 1998, Mardling & Lin 
2002) becomes 

Qp + g* (2) 
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For e < 1, 
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For the case in which only a synchronously spinning planet contributes to e, eqs (2) and (6) 
reduce to eq(l). 

Equations (2) and (3) allow us to examine the orbital evolution of star-planet systems 
with non-synchronized spins and arbitrary eccentricities. For example, a system in which 
the stellar tide can be neglected, (for instance, if the star were extremely compact), the tidal 
dissipation in the planet leads to eccentricity excitation if 



np>np,^§f,{e)n/Me)^§n 



(7) 



as e — > (also see Goldreich & Soter 1966). We now introduce two parameters which 
characterize the influence of the stellar tide and spin, as well as arbitrary eccentricity, on 
the orbital evolution of a system. The eccentricity damping efficiency factor {3 is defined as 
/5 = Tep/Te, where 

is the eccentricity damping timescale. For /3 < 1, stellar tides inhibit eccentricity damping, 
whereas for /? > 1 they enhance the damping rate. Negative values of {3 corresponds to 
eccentricity excitation and increasing the eccentricity always enhances the damping rate. 
From equations (1) to (3) we have 
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Here we have introduced a second parameter. A, to characterize the importance of the relative 
masses, radii and Q-values of the star and planet. Generally, tidal dissipation in the planet 
dominates the evolution when A is large. 



2.3. Magnitude of the Stellar and Planetary Q- values 

In eq (10), the value of Mp/M* can be directly determined from observational data (to 
within an uncertainty of sini, where i is the inclination of the system to the line of sight). 
The value of R*/Rp can be deduced theoretically from stellar and planetary structure models 
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or inferred empirically from the observed radius of the transiting planet HD 209458. The 
major uncertainty is associated with the magnitude of the Q-value because the dominant 
physical processes which determine the magnitude of * in both solar type stars and gas 
giant planets remain unresolved. Despite the availability of observationally inferred Q- values, 
it is not clear whether those obtained for solar-type binary stars and for Jupiter (see §1) can 
be directly applied to hot Jupiters. 

Theoretical analyses of tidal evolution mainly focus on the response and dissipation of 
both equilibrium and dynamical tides. The equilibrium model is based on the concept that 
a homogeneous spherical body continually adjusts to maintain a state of quasi-hydrostatic 
equilibrium in the varying gravitational potential of its orbital companion. Internal friction 
within the body induces the dissipation of energy and a phase lag which gives rise to a 
net tidal torque which transports angular momentum between the spin of the body and its 
companion and their orbits. Although turbulence can lead to dissipation in the extended 
convective envelopes of gaseous giant planets and low-mass stars, the convective turnover 
timescale is usually much longer than the period of the tidal forcing. Due to this reduction 
in the efficiency of the turbulent viscosity (Zahn 1977, 1989, Goodman & Oh 1999), the 
derived rate of angular momentum transfer falls short of that required for the circularization 
of solar-type binary stars by nearly two orders of magnitude (Terquem et al. 1998). For 
gas giant planets, turbulent dissipation of equilibrium tides gives Q ^ 5 x 10^^ (Goldreich & 
Nicholson 1977). 

In radiative stars, the tidal perturbation of the companion can induce the resonant exci- 
tation of low-frequency g-mode oscillations in the radiative region, which carry both energy 
and angular momentum fluxes (Cowling 1941; Zahn 1970). When this wavelike response is 
propagated to the stellar surface, radiative damping of the dynamical response provides an 
effective avenue for tidal dissipation in high-mass binary stars (Savonije & Papaloizou 1983, 
1984). However, the envelope of Jupiter is mostly convective (Guillot et al. 2003) and the 
relevance of g-mode oscillations is less well established (loannou & Lindzen 1993a,b, 1994). 
Nevertheless, the surface layers of hot Jupiters may attain a radiative state because they are 
intensely heated by their host stars. If so, g-mode oscillations may be excited just above the 
planet's convective envelope and dissipated through radiative or nonlinear damping (Lubow, 
Tout, & Livio 1997) as suggested for high- mass stars (Goldreich & Nicholson 1989). In rotat- 
ing gaseous planets, Coriolis force can provide a restoring force to the tidal perturbation and 
induce the resonant excitation of a rich spectrum of inertial modes (Savonije et al. 1995). 
When the forcing and initial-wave frequencies resonate with each other, the tidal response 
of the star is greatly enhanced (Papaloizou & Savonije 1997). Hot Jupiters probably formed 
as rapid rotators such that their tidal forcing occurs in the frequency range of inertial waves, 
which therefore constitute the natural and dominant response of the planet (Ogilvie & Lin 
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2004). 

In principle, the value of Q'^ ,^ depends in a highly erratic way on the resonant response 
which is a function of both the planet's spin and the perturber's forcing frequencies as well 
as the magnitude of viscosity. In the presence of viscosity, the spectrum of inertial waves is 
discrete, but as the viscosity tends to zero, the spectrum becomes increasingly dense (e.g., 
Dintrans & Ouyed 2001). In this limit, the dissipation rate may not simply vanish because 
of the increasing probability of resonant excitation of wave modes. Although the highly 
variable Qp values can introduce a dispersion in the rate of eccentricity dissipation, the 
spin of a planet changes as it undergo Kelvin Hclmholtz contraction such that the relevant 
frequency-averaged Q-value (~ 10^) may be asymptotically independent of the viscosity in 
the limit of small Ekman number (Ogilvie & Lin 2004). 

Based on these considerations and for computational convenience, we adopt the equilibrium- 
tide prescription in eq(2) for the determination of the eccentricity evolution of hot Jupiters. 



2.4. Eccentricity Evolution of Planets around Rapidly Spinning Stars 

The stellar tide can become stronger than the planetary tide if the spin rate of the star 
is high enough. From eqns(2) and (3) the critical stellar spin rate corresponding to e = is 



n-XQp, (11) 



Me). 

so that the eccentricity grows if fi* > f2*c- If flp ~ n and e -C 1, 

= (18 + 7A)n/ll. (12) 

For Jupiter-like planets orbiting solar-like stars, A = Xjq ~ 12 so that Q^c — 9n. 

In the next section, we show that planets attain a state of synchronization {Qp = n) 
on a timcscalc (tq^) which is generally much shorter than that of their host stars (r^^*). 
Wc therefore now consider systems in which the planet is nearly synchronously rotating. In 
Figure 2(a) we plot /? (eqn(9)) as a function of stellar spin period for various values of e. 
For illustration purposes we use A = Xjq with Qp ~ Q*. In addition, we have taken the 
orbital period to be 6.5 days, representing the three shortest period borderline planets. Now 
recall that /3 > 1 corresponds to enhanced eccentricity damping, < P < 1 corresponds 
to prolonged circularization times, while P < corresponds to eccentricity excitation. Also 
note that stars with radii the same as the Sun's and rotating with velocities = 50 km s^^ 
have spin periods of 1 day, while those with = 10 km s~^ have spin periods of 5 days 



-9- 



(see next section). For each value of fi*, /3 is bounded below by the value corresponding to 
e = (dotted curve). For high values of e, /? > 1 for most stellar spin rates, suggesting that 
eccentric borderline planets cannot have reached their present positions with eccentricities 
much higher than presently observed. 

Figure 2(b) plots (3 against fi* for the three shortest period borderhne systems HD - 
168746 (e = 0.08, M, = 0.92AfQ, AfpSini = 0.23Mj, orbital period P = 6.4 days, T^p = 
0.60-6.0 Gyr, A = 127.4; Pcpc ct al. 2002), HD 217107 (e = 0.14, M, = O.98M0, MpSini = 
1.30Mj, P = 7.1 days, r^p = 1.90 - 19.0 Gyr, A = 7.2; Fischer et al. 1999), and HD 68988 
(e = 0.15, A//, = 1.2Af0, Afpsin^ = 1.90A'/j, P = 6.3 days, r^p = 1.77 - 17.7 Gyr, A = 5.0; 
Vogt et al. 2002). The ranges for Tep are calculated using = Qp = 10^ — 10^ to represent 
the evolution of these parameters. The orbital evolution of HD 168746 was clearly rapid, 
since any eccentricity greater than its present value would enhance the damping rate even 
more (see Figure 2(b)). It is likely that the present value is, in fact, much closer to zero. 
Dissipation within the planet dominated the evolution because the planet-to-star mass ratio 
is relatively small, and the radii of Saturn-like planets are relatively large (Burrows et al. 
1997). Taking Rp = -Rsaturn = 0.84i?j (and sini = 1) gives the value for A quoted above 
In contrast, the eccentricity evolution of both HD 217107 and HD 68988 would have been 
retarded if the stars were rapid rotators initially. 



2.5. Dispersion in the Stellar Spin Frequency 

The spin frequency of a star is generally believed to be a function of its age, (see 
§3). However, in open clusters where member stars presumably have very similar ages, 
an order of magnitude dispersion has been observed in the stellar rotational velocity. For 
example, although most G dwarfs in the Pleiades are observed to have rotational speed 
VrSini — f2*i?^.sini < 10 km s"-*^, at least 20% of them are ultra-fast rotators (UFRs) with 
30 km s~^ < Ksini < 100 km s~^ (Soderblom et al. 1993a,b). The spin frequencies of these 
UFRs exceed the mean motion of planets with periods of several days. (For a star with a 
solar radius, for example, K = 30 km s~^ corresponding to a spin period 27r/Q* = 1.7 days.) 
In order to illustrate the importance of stellar spin, we carried out a Monte Carlo simulation 
in which a sample of planetary candidates was generated with the observed logarithmic 
period distribution between 3 days and 3 years. By adopting Mp = IMj, Rp = 1.35Rj, 
Af* = lAf0, R^ = IRq, Q'p = 10^, and Q'^ = 1.5 x 10^ and taking the planetary spins to 
be synchronous with the orbital motion, we deduced values for (3 and Tg given the observed 
spin periods of Pleiades stars and the simulated mean motions n — 2t:/ P. In Figure 3(a), we 
plot the inferred values of Te/r* (the age of the Pleiades cluster is taken to be 7 x 10^ yr) for 
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various simulated values of P. The open circles represent eccentricity damping timescales 
while the filled dots represent eccentricity excitation timescales. For all planets with P < 5 
days, the eccentricity is damped on timescales shorter than r^,. For borderline planets with 5 
days < P < 15 days, the eccentricities of several systems are damped on timescales Te > r^,. 
However, eccentricity is excited for several other systems in this period range. For P > 15 
days, the eccentricities of most systems are excited on timescales Tg ^ r^,. 

Ultra-fast rotators with V^sini up to 200 km s^^ have been observed among the stars 
of a Persei (Prosser 1992), an open cluster which is younger (r^, ~ 5 x 10 yr) than the 
Pleiades. Although the time available for eccentricity evolution is limited, the relatively 
small values (~ 1.5 x 10^ as inferred from observations of binary stars in these young 
clusters) together with medium to high stellar spin rates imply an enhanced role for the 
stellar tide and hence retarded eccentricity damping or even eccentricity excitation in some 
cases. Thus in these young clusters, we expect planets with periods greater than 1-2 days 
(including those with Tgp < r*) to have a noticeable dispersion in their period-eccentricity 
distribution. 

In older clusters such as the Hyades for which r* ~ 6 x 10^ yr, sin i is below 10 km s~^ 
for nearly all G and K dwarfs. In contrast, most F dwarfs are still fairly rapidly spinning, 
with rotation speeds in the range 20 — 70 km s~^. We performed a second Monte Carlo 
simulation, this time for the Hyades, adopting = 10^ and = 6 x 10^ yr. We found 
(Figure 3(b)) that the majority of the G and K dwarfs have slowed down sufficiently to 
promote eccentricity damping of short-period planets. There are nevertheless a few rapidly 
rotating F dwarfs which continue to promote eccentricity excitation. 



2.6. Stellar and Planeteiry Structural Adjustments 

Eqn (2) indicates that e is a rapidly increasing function of Rp. The evolution of a and 
Qp^^ also depend sensitively on Rp and R^, (see the next section). During the early stages 
of stellar evolution, both stellar and planetary radii are somewhat larger than their values 
at later stages. We now assess the effects of stellar and planetary structural adjustment in 
regulating eccentricity and semi-major axis evolution. For computational convenience, we 
approximate existing numerical stellar models (Iben 1965, 1967) with the analytic formula 

R, ~ P0(r,/lOV)"°-' (13) 

for T* < 10^ yr. This timescale for solar-type stars to undergo Kelvin Helmholtz contraction 
onto the main sequence is much shorter than Tg or tq*, such that the initial large stellar radii 
cannot significantly modify e or Qp. 
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However, planets have much slower Kelvin- Helmholtz contraction rates, (Burrows et al. 
2000, Bodenheimer et al. 2000), and are such that 



Rp ~ 1.2i?j(r,/ro)-°-\ (14) 

where Tg — 1 Gyr. For r* = To, Rp in this prescription is actually smaller than that used 

in our models above {Rp = 1.35Rj). The larger value takes into account stellar insolation 
of a short-period planet's surface, a process which tends to extend its contraction timescale 
(Burrows et al. 2000) and tidal inflation (Bodenheimer et al. 2001). Using eqs(2) and (14) 
for an unheated planet, we flnd that is an order of magnitude smaller when = 10® yr 
than when T^, = 4.6 Gyr. (For = IQ^-^^*^-^ yr, a planet's size is essentially independent 
of its mass.) Nevertheless, during the first 10® yr, Te/r* > 1 and there is insufficient time 
to circularize the orbits of borderline planets. However, since Tep/r* oc t~°'^, eccentricity 
damping occurs mostly during the later stages of main sequence evolution and borderline 
planets can be circularized provided their host stars are not rapid rotators. 

In addition to Kelvin Helmholtz contraction, a planet's size is modified by tidal dissi- 
pation which leads to an internal energy generation rate 

Et = + apMpRlnppp] (15) 

where \e\p refers to the planet's contribution to the circularization process in eq(2) and Qp is 
the rate of change of Op (see next section). For sufficiently large \e\p or \Qp\, the magnitude 
of Et exceeds the rate of energy released ExniMp, Rp) due to an unheated planet's Kelvin- 
Helmholtz contraction. In this limit, a planet's radius evolves toward an equilibrium value 
Re such that its surface luminosity, C{Mp, Re), is balanced by Et (Bodenheimer et al. 2001). 

However, the evolution of Rp in eq(14) depends on Ekh- If a short-period planet arrives 
at its current location shortly after formation, the value of E^h at the time will exceed Et. 
After that, Kelvin-Helmholtz contraction of the planet will stall once E^h becomes less than 
Et. On the other hand, if a short-period planet acquires its present orbital properties after 
it contracts to a radius Rp < Re, its size will adjust at a rate (Mardhng & Lin 2002) 

rtnUit 



with an enlargement timescale 



Rp M^Rp . . 
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This timescale is generally longer than TQp, in which case a planet cannot inflate significantly 
before it attains a state of near synchronization (Gu et al. 2003). 

But, interior to a ~ 0.04AU, tr < T^p for a planet with non-negligible eccentricity 
(> 0.1) so that it can infiate before the eccentricity is damped out. Planetary infiation 
enhances the rate of tidal dissipation and shortens Tg. For some values of e and a, Re 

exceeds a planet's Roche radius, Rr — (Mp/3M*)^/"'a, so that it begins to lose mass through 
Roche lobe overflow while its radius is confined to Rp = R^ (Gu et al. 2003). Exterior 
to 0.04 AU, the condition required for planetary infiation {tr < r^p) is more difficult to 
satisfy unless the actual eccentricity- damping timescale, Tg, is significantly prolonged by 
eccentricity-excitation due to the tidal dissipation inside a rapidly rotating host star, or due 
to the secular perturbation of other planets. Nevertheless, is a rapidly decreasing function 
of a. Thus for a borderline planet, if a is sufficiently large and Et sufficiently small, R^ is 
unlikely to be much larger than our adopted values of 1.35 — l.bRj. 

2.7. Eccentricity Dispersion Among Borderline Planets 

Among the known planets, the shortest period planet with a significant eccentricity is 
HD 68988, with e = 0.15, P — 6.276 days, and minimum mass MpSini — 1.9Mj. The mass 
and age of the host star are = 1.2 Mq and r* ~ 6 Gyr respectively. Given these values 
and taking Rp = 1.35i?j, T^p ~ n if Q'p ^ 10^ (Vogt et al. 2002). The value for Rp is that 
inferred from transit observations of the short-period planet around HD 209458 (Brown et 
al. 2001). 

The dynamical properties of HD 68988 are remarkably similar to those of HD 217107 
which has a similar eccentricity (e = 0.14). However, there are several other planets such 
as HD 168746 and HD 130322 which have similar values for Tgp but negligible eccentricities. 
We propose that the dispersion in eccentricity of borderline planets is correlated with the 
dispersion of stellar spin frequencies of young stars. Based on the lack of chromospheric 
activity, Vogt et al. (2002) inferred that HD 68988 is a slow rotator. If correct, the stellar 
tide will be promoting eccentricity damping in this system today. Using the Skumanich law 
(see eqn (34) below) as well as data from Soderblom et al. (1993a, b), we find that the mean 
spin frequency of F, G, and K stars with an age comparable to that of the Pleiades cluster is 
about half the present orbital mean motion of the planet around HD 68988, i.e. Jl* < Jlpc- 
However, as discussed above, slow and fast rotators coexist in this cluster, with the dispersion 
in the observed values of Vr sin i being large. There also exist solar-type stars which are older 
than the Sun but which are more rapidly spinning. If HD 68988 and HD 217107 were rapid 
rotators when they were young, at that time /3 < 1 so that Tg > T^p (Figure 2(b)). The effect 
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of a young star's tide is also enhanced by its slightly smaller value of Q'^ as inferred from 
observation. Similarly, if HD 168746 and HD 130322 were slow rotators with < since 
birth, Te < Tep. 

To support our conjecture, we carried out some numerical simulations using a scheme 
which calculates the tidal, spin and dynamical evolution of a multi-planetary system (Mard- 
ling & Lin, 2002). Using the observed kinematic data for HD 68988 together with Hp — 
l.5Rj. = 1.3i?0, Q'^ ~ ~ 10^, = n, and fi* = 27r/(2 days), we obtain Tg ~ r^, ~ 
a few Gyr. However, if the spin frequency of the star is halved, Tg is reduced by a similar 
factor. The magnitude of becomes less than 1 Gyr for slowly rotating host stars with 
< n. These numerical results are in general agreement with our conjecture. A more 
comprehensive numerical analysis is described at the end of the next section. 



3. Evolution of the Stellar and Planetary Spins 



The above results clearly demonstrate the importance of the stellar tide and spin in 
the eccentricity evolution of short-period planets. It seems reasonable, therefore, that the 
dispersion in the e — Tep relation is a result of the spread of spins among young stars. 
However, most of the target stars in the radial velocity surveys are chosen for their low 
level chromospheric activities. These type of stars, including several borderline systems with 
measurable eccentricities, have slow rotation speeds like the Sun. 

Nevertheless, some, if not all, host stars of extrasolar planets may have been rapidly 
spinning in the past. It is therefore worthwhile investigating the evolution of Jl* and ftp. 
Assuming an aligned rotator and that ^ Mp, the rate of change of the stellar spin 
frequency is given by (Mardhng & Lin 2002) 

,2 \ /7i.r\2/n\3 



n 



Mp 



/3(e) -/4(e) 



n 



(18) 



where 



/3(e) = (l + fe^ + fe^ + Ae«)/(l-eY, (19) 
Me) = (l + 3e^ + |e4)/(l-eY/^ (20) 

and UJ^, is the rate of change of the stellar spin due to angular momentum loss via a stellar 
wind (see the next subsection). The quantity e* is the stellar mass fraction participating in 
tidally induced and external angular momentum exchange and a* is the moment of inertia 
coefficient of the star. The rate of change of the planet's spin is given by 

,2 \/,^\/n\3 



n 



-pOipQp 



Rr 



Me) - Me) 
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+ Co. 



(21) 
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where and ap are the planetary counterparts of e* and a*. Note that depends more 
strongly on the mass ratio than does Clp. Since gaseous planets have extensive convection 
zones which can be fully mixed, ~ 1. In contrast, solar-type stars have shallow surface 
convection zones so that is a few times 10~^ for G dwarfs, and an order of magnitude 
smaller /larger for F/K dwarfs. 



3.1. Timescales for Spin and Semi-Major Axis Evolution 

The expression for the rate of change of the stellar spin (eq (18) ) involves two contri- 
butions, as does its counterpart for the rate of change of planetary spin. For post-formation 
evolution we can neglect ujp for planets, while cj* is associated with angular momentum 
loss via a stellar wind. The tidal transfer of angular momentum induces Vlp to evolve on a 
timescale 

for non-synchronously spinning planets (see eqns (1) and (18) with e ^ 1 and M^, ^ 
Mp). Therefore the planetary synchronization process is much more rapid than the orbital 
circularization process. However, as a consequence of this, the planetary spin can only be 
guasi-synchronized with the orbit while the eccentricity is non-zero. From eqn (21), this 
value is Q,p = h{e)n/ fi{e) — > (1 -I- 6e^)n for e ^ 1. 

For conservative systems {uj^, = 0), the stellar spin evolves on a timescale 

" = r-''^ (t) (t) 

so that for most systems of interest (in the absence of angular momentum loss) we have 
that tq^ > Tep. Thus we can conclude that fi* does not change significantly prior to the 
circularization of a planet's orbit (in the context of solid satellites around gaseous giant 
planets, the relative values of tq^^, t^, and tq* are even more extreme because Q' values 
for satellites arc much smaller than those for such planets; Goldrcich & Soter 1966). We 
conclude, therefore, that given the slow stellar rotation rates observed today in borderline 
systems, were they to have suffered no angular momentum loss in the past, tidal dissipation 
inside the stars could only ever have enhanced eccentricity damping. 

Assuming that the spins of a planet and a star are aligned with the orbit normal of such 
a system, the total angular momentum perpendicular to the orbit is 



-'total (a, fi*, e) = Jo{a, e) + Jp{a) + J*(Q*), 



(24) 
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where Jo — MpO^ny/l — is the orbital angular momentum. The angular momenta of the 
planet and star are 

Jp^apMpRlQp (25) 

and 

J* = a^e^M^Rl (26) 

respectively, where we have assumed that the whole planet participates dynamically whereas 
only a fraction e* of the star is involved in angular momentum exchange. The rate of change 
of the total angular momentum is 

— - ^— ^ j + apMpRl Qp + a^e^M^Rl (27) 

where 

j^, = a^e^M^Rlu^ (28) 

is the rate at which the system loses angular momentum via a stellar wind. In a conserved 
system, we find from Eqns (18) and (21), 

to first order in the eccentricity. From Eqns (1) and (8), this can be written as 

d /2e2 



where 

in the limit that i7p = n. For small to moderate eccentricities and for most stellar spin periods 
of interest, as well as for a wide range of star-planet masses and radii (characterized by A), 
~ 0(1) (Figure 2). However, as illustrated by the shortest period borderline planets 
described in Section 2.2, A can vary widely being ~ 0(1) for systems with relatively massive 
planets and lower-mass stars, to several hundreds for systems with lower mass planets and 
higher mass stars. For systems in which f2,, > n, d > as long as < a (5/2. On the other 
hand, d < whenever < n. 

For Jupiter-mass planets around solar-type stars, \Jp\ ^ |J*| so that the stellar spin 
provides the dominant source of angular momentum exchange with the planet's orbit. In 
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the absence of a stellar wind, stable configurations satisfy Jo > 3 J* (Hut 1980, Lin 1981), or 
equivalently. 



a > ar 



3Q;*e* 



n 



1/2 



(32) 



y(l-e2) V^P. 

and angular momentum exchange acts to reduce the difference between n and until an 
equilibrium configuration is achieved. The timescale for the semi-major axis to evolve is 

2 



a 

\a\ 



In - 



(33) 



3.2. Simultaneous Evolution of the Stellar Spin and Planeteiry Orbit 

Stars are known to lose spin angular momentum via stellar winds so that ci;* in Eqn (18) 
is finite and negative. The youngest stars are observed to have rapid rotational velocities 
which are sometimes as high as their breakup speeds (Stassun et al. 1999). Although 
observations of G and K dwarfs in the a Persei and Pleiades clusters (with r* < 10^ yr) 
indicate a large dispersion in Ksini (with magnitudes ranging from a few to 100-200 km 
s~^), the same type of stars in the Hyades (with r* = 6 x 10^ yr) have nearly uniform Vr sini 
which are less than 10 km s^^ (Soderblom et al. 1993a,b). The simplest description of the 
decline in the mean rotation velocity of mature G and K dwarfs is given by Skumanich (1972) 
such that 

< K sini >~ V;(r*/ro)-°-^ (34) 
where ~ 4 km s~^ and Tg = 1 Gyr. A more general description is given by 



Vo Y f Vo 



V^(t*) sin i ) \ ^1 sin i 



(35) 



which more accurately models both the dispersion in among young stars at some initial 
time Ti and the homogeneous low K among mature stars with r* ^ t\ and To (also see 
Mayor & Mermilliod 1991 who use a different power index). For G and K dwarfs, t\ ~ 10^ 
yr and V\^m.i — V^(ri) sini ~ 30 km s~^. However, the F dwarfs in the Hyades lose very 
little angular momentum through stellar winds and sin i continues to extend up to 100 
km s~^ at r* = 6 X 10^ yr. Thus Vq, l^i. To, ti, and the power indices in eqs (34) and (35) 
may be functions of M*. 

From the empirical equation (34), we deduce that 
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where differentiation is with respect to r*, = and 7 is a calibration factor which is 

set to unity for G and K dwarfs. Eq(35) yields a similar expression: 



2 N -1-5 



2R,To\\V,smiJ To J ^ ^ 

In both expressions, the timescale for spin-down induced by a stellar wind is 

T,, = Vr/\Vr\ ^ {2Toh){no/^*f. (38) 



Due to the initial rapid rate of angular momentum loss, r^, <^ r^^ for rapidly rotating 
young G and K dwarfs with greater than a few times 10^ yr (i.e the wind-loss timescale is 
much shorter then the tidally driven synchronization timescales). The evolution of the spin 
rates of these stars is mainly determined by eq(34) rather than eq(18). However, after these 
have declined to sufficiently small values, the establishment of a quasi-equilibrium state is 
possible in which the rate of angular momentum loss via a stellar wind (eq 28) is balanced 
by the rate at which the star gains angular momentum from the orbit as the planet attempts 
to spin up the star. Under such circumstances, = so that from eqn (18) an equilibrium 
spin rate 

^e = #|n. + ^^^^W7M^^ (39) 



/4(e) 9 \U{e)nJ \MJ \R 
is possible provided that 

|.„.I<H«£LA. (40) 

7 VI - e"' 

In eqn (39), Oe is the equilibrium semi-major axis and rie — n{ae) is the associated mean 

motion. This quasi-equilibrium can be established through the evolution of either or a 
so that it is attainable on a timescale given by the minimum of r^^, , tq^^, and r^. For a given 

total angular momentum, ii flp = n then the equilibrium quantities of Og and can be 
determined self-consistently from eqs (39) and (24) as functions of e. 

Jtotai(ae, CO., e) = al/^A{e, cu,) + al'''B{e) + af'Cie) (41) 

Aie.uj^) — — — , , (42 

^ ' ' 9U{e)VGM:M^R,' ^ ' 

B{e) = Mp{GM,{l-e')y'\ (43) 

C(e) = apMpRly/GAf, + a,e,M,Rlh{e)^/GK/ U{e). (44) 



where 



and 
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In Figure 4(a), we plot ag and as functions of Jtotai for three values (0, 0.2, 0.5) of 
the eccentricity. The total angular momentum is normalized by Jq.iau, the value of Jupiter's 
angular momentum in a O.IAU orbit around the sun. Similar to the calculation presented 
in Figure 2, we adopt M, = IMq, = IRq, Mp = IMj, Rp = 1.35i?j, Q'p = Q'^ = 10^ 
a^,e^ = 0.01, and apCp = 0.2. However, the establishment of this equilibrium for a wide range 
in a is only possible for a very slow rate of spin-down, so we set ci;* to be 10""^ that in eq(36). 
In Figure 4(b), we represent a slower spin down rate with the same model parameters except 
we set a;* to be one tenth that in Figure 4(a). 

The results in Figures 4(a)and 4(b) show that is sensitive to both e and Jtotai, while 
fie varies primarily with Jtotai- The structures of the curves are indicative of the dominant 
source of angular momentum in the system, and thus the behavior of the equilibrium is highly 
sensitive to the planets position. For low eccentricity systems with > 0.1, the dominant 
source of angular momentum is the orbit. Because Jtotai — Jw < 0, angular momentum is 
continually drained from the system, but this only requires minute adjustments to Og. For 
low eccentricity systems with Og < O.lAU the stellar spin is more prominent and in order 
to maintain this spin, larger changes in Og are necessary. When ag < O.IAU, a planet will 
move closer to the star as its orbit is circularized, and this will be accompanied by the 
spin-up of the star. Based on these results (with the planetary and stellar masses and radii 
given earlier) we infer that as angular momentum is continually lost via stellar winds, the 
state of near synchronization (in which Qp = n ~ Jl* and e ~ 0) occurs at a semi-major 
axis of approximately 0.02-0.06 AU and Jtotai/ Jo.iAu ~ 1-0, corresponding to a total angular 
momentum of approximately lO^^g cro? /s. It is approached asymptotically as Jtotai slowly 
decreases. 



3.3. Present Spins of Host Stars 

In general there are two classes of quasi-angular momentum equilibria. The close equi- 
libria with Cg < Cg are unstable to tidally induced angular momentum exchange because they 
enhance the difference between and n and lead to further orbital evolution (Hut 1980, 
Lin 1981, Rasio et al. 1996). There are also those distant equilibria with ag > Qc which are 
stable against tidal interaction because they lead to the reduction of difference between 
and n. 

The magnitude of Og does not explicitly depend on Jtotai (see eq 32). However, when 
the stellar spin equilibrium is attained, its magnitude ago is determined by and n, which 
along with Og and Jtotai are functions of e. For the above model parameters, Qco = 0.029AU 
for e = 0.2. We plot the value of Ogo and the corresponding Jtotai for various values of e in 
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Figures 4(a)and 4(b). 

The results in these figures also show that is a decreasing function of J total (eq 41). 
In the large Jtotai limit, if most of the angular momentum is retained by the planet's orbit, 
a would be large and the rate of tidal transfer of angular momentum would generally be too 
small to balance the spin down induced by the stellar wind. However, if a large amount of 
angular momentum is stored in the spin of the planets and their host stars, a small a^. would 
be sufficient to balance against the large magnitude of J^^. For modest values of JtotaU it is 
possible for angular momentum to be mostly contained in the planet's orbit and for the host 
star to spin relatively slowly with a small magnitude of J^. In this case, a spin equilibrium 
can be attained with modest values of Ue- Note that there is an upper limit to the magnitude 
of ae (~ 0.2 — 0.3AU), larger than which the rate of tidal transfer of angular momentum 
cannot be balanced with the wind loss rate. 

The evolutionary fate of planets depends on the magnitude of Jtotai when a reaches Og. 
Consider the possibility that a planet has migrated, through its interaction with a nascent 
disk of gas, residual planetesimals, or other planets, close to the surface of a rapidly spinning 
host star. If the planet is deposited in a region where Jl* is smaller then n the planet will 
not survive. Transfer of angular momentum from the orbit to the stellar spin will cause 
the planet to migrate into the star. However, if the planet is deposited outside co-rotation 

> n) its orbit will expand through star-planet tidal interaction while the star spins down 
through wind losses until the system comes into equilibrium at a = a^.. As the star continues 
to loose angular momentum through its winds, the total angular momentum of the system 
will decrease and the location of Og will move inward, as can be seen in Figures 4(a) and 
4(b). If the system still retains a relatively large Jtotai (stored primarily in the stellar spin) 
the magnitude angular momentum lost through winds will be too great for the system to 
remain in equilibrium. The planet's inward migration would continue until it merges with 
its host star on a time scale Ta(ae) < Tn*(ae) which, for planets with P < 2 days, is less 
than a few Gyr. However, if a spin equilibrium is retained after Jtotai has already reduced 
to modest values such that > ^co, it would be stablized against further evolution. Note 
that the values of a^o correspond to orbital periods in the range of 3-4 days, which are close 
to the observed cut off in the period distribution of extra solar planets. 

In the previous subsection, we indicated that main sequence stars continually lose angu- 
lar momentum through stellar winds. For example, the ultra fast rotators among the young 
G and K dwarfs spin down rapidly. During their spin down, the of these stars may become 
comparable to the n of their short-period planets. But this state of near synchronization 
cannot be maintained because the rate of tidal transfer of angular momentum between a 
star and its planet, J*, is inadequate to balance that of wind-driven loss of the stellar spin 
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angular momentum J^j. Consequently, reduces to D,e << n and a declines through a 
sequence of non-synchronized (with < ~ n) quasi-equilibria on the timescale of To. 

This evolutionary path may have been taken by some G and K dwarf host stars of the 
shortest-known-period (~ 3 days) planets such as HD 46375 (Marcy et al. 2001), HD 177949 
(Tinney et al. 2001), and HD 187123 (Butler et al. 1999). The corresponding K for 
synchronous rotation (with Q^, = n) in these systems is ~ 15 km s~^ which is comparable to 
that of typical dwarf stars in the a Persei and Pleaides clusters (assuming the expectation 
value < sini >= 0.7). It is entirely possible that some time in the past, the spins of these 
stars may have been synchronous with the mean motion of their planets. With their present 
configurations and Q'^ = 10^, their present r^* of a few Gyr is comparable to their estimated 
T*. Yet the spin of these stars is slow (Q* < n), as indicated by their quite chromosphere. 
We suggest that the effect of the tidal interaction is overshadowed by that of the stellar wind 
in these G and K dwarfs when they evolved through a state of synchronization. Although 
they have attained a stable quasi-equilibrium today with a ~ 2ac, their orbit may decay into 
an unstable configuration with a < Qc- In order for these planets to have survived ~ a 
few Gyr, their Ta ~ 2rn* > which requires Q'^ to be comparable to or greater than 10^. 
In principle, a nearly synchronized quasi-equilibrium (with ~ ~ n) is attainable for 
the longer-period (with P > 10 days) planets around G and K dwarfs, but the timescale for 
evolution into such a state is much longer than the age of the mature stars {i.e. rn* >> r*) 
which makes their potential existence coincidental. 

But for the more massive F dwarfs, r^* < r^^* because their u;* is much smaller than 
that of the G and K dwarfs. For these stars, we can use eq (36) for a;^. by setting 7 ~ 0.1. In 
this case, the state of quasi-equilibrium may be maintained because the transitory decline 
of Jl* can be compensated by a tidal transfer of angular momentum from the planet's orbit 
to the stellar spin, resulting in an net increase in both and n{ae) and a reduction in Ue- 
In addition, the relatively shallow convection zone of the F dwarfs also reduces their and 
Qc well below that of the G and K dwarfs. It is thus possible to attain a sequence of nearly 
synchronized (with ^2* ~ ~ n) quasi equilibria with > a^. This sequence of events may 
have led to the present orbital configuration of the planet around r Boo. 

3.4. Numerical Models of Spin-Orbit Evolution 

In order to verify the above qualitative and analytic description, we numerically inte- 
grate, with a fourth-order Runge-Kutter scheme, eqs (2), (18), (21), and (30) with appropri- 
ate substitutions for /i, 2,3,4, 1^*, Jp,*- The dependent variables are flp^^,a, and e. We adopt 
the standard model parameters in which = IMq, Mp = IMj, = IRq, Rp = l.5Rj, 
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Up = 0, and Q'j, = 10^. Assuming that Q'^ is around 1.5 x 10^ before t = 10^ yr and 10^ after 
t = 10^ yr, the star's dissipation parameter is taken to be of the form Q'^ — a+6erf [(r'— 5)/3], 
where r' = yr, a = 5.75 x 10^ and 6 = 4.2 x 10^ 

In order to demonstrate several possible outcomes, we constructed six models. In models 
1 and 2, we set P = 3 days at i = to represent the shortest period planets. In models 3 
and 4, P is set to be 7 days at t = as a possible initial condition of borderline planets. In 
models 5 and 6, we choose P = 1.5 days to represent some hypothetical ultra-short-period 
planets which may have once existed close to the surface of their host stars. To explore the 
effects of stellar spin, models 1, 3, and 5 start with = lOn, while in models 2, 4, and 6, 
f2* = 9n. These spin rates are rapid, but not near the break up values. In all models, we 
choose Qp = lOn and e = 0.2 at i = 0. All models are integrated for 6 Gyr. 

In Figure 5(a), we plot the e evolution of models 1, 2, 3, and 4 which are represented 
with solid, short-dashed, long-dashed, and dot-dashed lines respectively. Both models 5 and 
6 are represented by dotted hues. In Figures 5(b), 5(c), and 5(d), we plot the evolution of 
a, flp/n and fl^/n respectively. 

After brief initial episodes of e excitation, e declines rapidly in models 1, 2, 5, and 6. 
Close to the star, Tg ~ 3 Gyr for the short-period planets, and rapid late stage e damping 
erases any trace of early e excitation. Despite an initial decline in the eccentricities in 
models 5 and 6, the planet migrates outward and e levels off. Eventually it is caught at in 
the quasi-equilibrium at a^, and e begins to drop again. In models 3 and 4, the eccentricity 
is maintained at non-zero values even after 6 Gyr. The contrasts seen in the models with 
different stellar spin rates supports our conjecture that the dispersion in the e—P distribution 
may be due to tidal dissipation in host stars with various rotation speeds. 

The planets quickly decrease in spin frequency due to the artificially high initial values. 
They are all able to attain synchronous spin frequency within ~ 1 — 2 Gyr, with the closer 
planets doing so more quickly. The effect of the stellar spin on planetary spin is negligible. 

After 400 — 500 Myr, all the models are caught in the quasi-equihbrium described earlier 
in this paper. As angular momentum is lost to stellar-winds it is balanced by tidal transfer 
from the orbit to the star. Both the semi- major axis and stellar spin reach their equilibrium 
values described by eq (39), eq (41) and shown in Figure 4. Models 2, 4, and 6, that start 
with slower initial stellar spins, have correspondingly smaller equilibrium stellar spins and 
semi-major axes. Any subsequent evolution proceeds on much longer timescales governed 
by LJ*. 

In all the models the stellar spin steadily decreases from a large initial value. However, 
Figure 5(d) reflects changes in both stellar spin and semi-major axis. The stellar spin in 
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model 1 initially decreases more quickly then the orbital frequency, but by ~ 100 Myr the 

orbital frequency is decreasing faster and Vl^/n increases. The same behavior occurs in model 
5, but at an earlier time. Models 3 and 4 have such large semi-major axes that the stellar spin 
must decrease significantly before any equilibria can be established. In Models 3, 4, 5, and 6 
the stellar spin has stopped changing by the time the semi-major axis reaches its equilibrium 
value. In contrast, the stellar spin in models 1 and 2 continues to decrease for about 100 
Myr, even after the semi-major axis has stopped changing significantly. Further evolution 
should slowly move the system from quasi-equilibrium into synchronization (Q^ = n ~ 
and e ~ 0). 



4. Summciry and Discussion 

In this paper we have considered the effect of tidal interaction between main sequence 
dwarf stars and their short-period planets. It is generally customary to assume that the or- 
bital eccentricity of planets is primarily damped by the dissipation of stellar tidal disturbance 
inside the planets. Such a scenario provides a reasonable account for the circular orbits of 
planets with periods less than about a week. But for planets with periods slightly longer than 
a week, there is a spread of orbital eccentricity. We suggest that this dispersion is associated 
with different past rotational speed and a range of spin down rate among their host stars. 
This conjecture can be best verified by the detection of short-period planets around stars in 
young clusters such as a Persei and Pleiades. We expect short-period planets around rapidly 
rotating F, G, and K dwarfs to have relatively large eccentricities. For older clusters such 
as the Hyades, short-period planets around F dwarfs should have larger eccentricities than 
those around G and K dwarfs, mainly due to the difference in the stellar spin-down rate. 

Rapid stellar spin also causes outward orbital migration on slightly longer timescales. 
Thus we expect the minimum orbital period for planets around rapidly rotating dwarfs in 
young clusters to be less than its present 3 day value. In contrast, this cutoff period for slowly 
rotating dwarfs in young clusters should be the same as its present value. It is also possible 
that the period cutoff represents the boundary corresponding to survivability against tidally 
induced orbital decay. In that case, the period cutoff for slowly rotating dwarfs in young 
clusters would be shorter than 3 days. 

During their spin down, the stars' spin angular frequency temporarily coincides with 
the mean motion of the planets. Among the short-period systems, this state of near syn- 
chronous rotation can only be preserved around F dwarfs which spin down gradually. All 
main sequence dwarfs attain a quasi-angular momentum equilibrium in which the host stars' 
loss of angular momentum through stellar winds is balanced by their tidal transfer of angular 
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momentum with their planets. For G and K dwarfs, the equihbrium frequencies are much 
smaller than that for F dwarfs. Consequently, it is much easier to form short-period nearly 
synchronous systems around F dwarf than around G and K dwarfs. 

Despite the enormous progress that has been made in the area of extra solar planet 
searches, there is insufficient data to verify out conjecture. But future radial-velocity, astro- 
metric, and transit searches for planets in clusters of various ages should provide sufficient 
data to test the scenario presented here. 
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Fig. 1. — (a) Period-eccentricity relation for extrasolar planets, (b) Eccentricity as a function of 
eccentricity-damping timescale for planets with orbital periods less than two months, (c) Eccentric- 
ity as a function of eccentricity-damping timescale normalized by the stellar age, (d) Eccentricity 
as a function of stellar flux received at planets orbital radius. 
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Fig. 2. — Eccentricity damping efficiency factor /3 vs. stellar spin period, (a) Sun/Jupiter param- 
eters with Q* = Qp, P = 6.5 days and various eccentricities, (b) Three representative borderline 
systems. 
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Fig. 3. — Te/r^ for stars as a function of orbital period. Empty circles correspond to eccentricity 

damping and filled circles are excitation timcscalcs. (a) Pleiades cluster (r* = 70Myr). (b) Hyades 
cluster (r^ = 600Myr). Tg was calculated assuming planet-orbit synchronization and eccentricities 
between 0.0 and 0.7 (see text for more details). 
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Fig. 4. — (a) The equilibrium semi-major axis and stellar spin period as a function of total angular 

momentum in the system for three eccentricities: e=0 (solid), c=0.2 (dotted), and c=0.5 (dashed). 
The equilibrium is reached when the angular momentum lost to stellar winds is balanced by that 
transfcred from the orbit and spin of the planet such that the stellar spin remains constant, (b) 
Stellar wind set at one-tenth that in (a). 
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Fig. 5. — Results of a long-term integration of eccentricity, semi-major axis, stellar spin frequency 
and planetary spin frequency. Models 1, 2, 3 , and 4 are represented with solid, short-dashed, long- 
dashed, and dot-dashed lines respectively. Models 5 and 6 are represented with dotted lines, (a) 
Eccentricity as a function of time, (b) Semi-major axis as a function of time, (c) Ratio of planetary 
spin to orbital frequency as a function of time, (d) Ratio of stellar spin to orbital frequency as a 
function of time. 



